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Abstract 



In this paper we derive a formula for the density of states in the presence 
of inelastic scattering in the quantum well of a double barrier structure as 
a function of a characteristic time of the motion of electrons (namely, the 
round trip time in the well) and of transmission probabilities for the whole 
structure and for each barrier. In the model we use the scattering processes 
due to phonons, impurities, and interface roughness, are taken into account 
by a unique phenomenological parameter, the mean free path, which plays the 
role of a relaxation length. We also show that, for lower rates of incoherent 
processes, the derived formula reduces to the one obtained by means of the 
Breit-Wigner formalism. 

PACS numbers: 71.20.-b, 73.20.Dx, 73.40.Gk 
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Typeset using REVT^ 



I. INTRODUCTION 



The density of states is one of the most important quantities for the study of equihbrium 
and transport properties of quantum effect devices. A recently derived formulaS estabhshes 
a simple general relation between the density of states in a mesoscopic system and the dwell 
times for each incoming channel connecting the system to the external world. 

In this paper we apply that result to the quantum well of a double barrier resonant 
tunneling diode. Moreover, in order to get closer to real systems and experimental results, 
we consider the effects of inelastic processes taking place in the well, which are not accounted 
for in the general relation of Ref. |I[ 

We obtain a very compact formula connecting the density of states in the well to a 
characteristic time of electron motion in double barrier structures, i.e., the time an electron 
takes to complete a round trip of the well, and to the transmission probabilities for the 
whole structure and for each barrier. 

In the literature the density of states in the well is usually obtained using Breit-Wigner 
formulas,i^ which lead to very simple and compact expressions. Dissipation can be easily 
accounted for by introducing a partial resonance width for all the inelastic processes.ui 
However, for Breit-Wigner formulas to hold true, it is necessary that all the partial resonance 
widths be much smaller than the separation between the resonant energy levels, and between 
each level and the top and bottom of the potential well. This condition establishes an upper 
limit on the rate of incoherent processes for the applicability of the Breit-Wigner formulas. 

On the contrary, our model is valid even when coherence is completely destroyed. Scat- 
tering with phonons, impurities, and interface roughness is accounted for by a unique phe- 
nomeno logical parameter /, the mean free path, a concept which is well established in solid 
state physics.0 

We assume that an electron traversing an infinitesimal length dx of the one-dimensional 
device structure, experiences a collision with probability dx/l, and that electrons emerge 
from collisions with an equilibrium distribution function in a state with completely random 
phase. This means, for instance, that the square amplitude of a plane wave function of wave 
vector k attenuates exponentially as it propagates along the x-axis, with a characteristic 
length equal to /. It is a situation very similar to an electromagnetic wave propagating in 
a dissipative medium. The difference is that the number of electrons has to be conserved, 
therefore electrons which seem to have disappeared have actually made a transition to a 
different state, with a phase completely uncorrelated, so that there is no quantum inter- 
ference between these electrons and those that have not undergone an incoherent process. 
In this model, all collision processes are effective in randomizing phase and energy, and we 
do not make any difference between the effects of elastic scattering (due to impurities and 
interface roughness) and inelastic scattering (with phonons, for instance). A more sophisti- 
cated model should take into account these differences, and, as a minor improvement, phase 
randomization and energy relaxation could be split using a different characteristic length 
for each process. 

BiittikeiS proposed a model for the inclusion of incoherent processes which is similar 
to the one we use. There is an inelastic scatterer in the well modeled by an extra branch 
leading away from the conductor to an extra reservoir, which does not draw net current, but 
permits phase randomizing events. Anyway, such a model is valid for very small differences 
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between electrode chemical Dotentials and/or when energy relaxation is not accounted for.S 
In the model by Knabchenp the inelastic scattering probability e for an electron traversing 
the well introduced by Biittiker is simply substituted by exp(— w//), where w is the well 
width and / the mean free path. 

In our model scattering is spread over the whole region, and not concentrated in a single 
point, and any potential profile can be considered. No addictional condition is required to 
obtain formula for the density of states. The hypothesis of smooth potential in the 
well is imposed in order to obtain the compact formula (0). The importance of the energy 
relaxation mechanism will be shown elsewhere.0 

As we shall show in section III, we use in our formula a characteristic time of electron 
motion in the well defined on the basis of Larmor times for transmission and reflection.EUll 
The tunneling time problem is the subject of a long-standing controversy: while the dwell 
time0 is widely accepted in the scientific community, there is no consensus on the actual 
time spent in a region by transmitted and reflected particles,liiliii2l due to both the fact 
that there is no operator for time in quantum mechanics (therefore we cannot perform 
a direct measurement of time) and that electrons do not follow actual trajectories in the 
Copenhagen interpretation.!! The Larmor times are obtained as the result of an indirect 
measurement: a weak perturbation is applied to the region of interest (i.e., a magnetic field, a 
real potential, or an imaginary potential) and some variation in the properties of transmitted 
and reflected particles is measured (spin precession, phase rotation, or particle absorption, 
respectively) .llHll'0 What is controversial about Larmor times is the interpretation of such 
results of an indirect measurement as the "actual" times spent in the considered region. 
However, this point is not relevant to the aim of the present work, where we are just interested 
in deriving a relation between the Larmor times and the density of states in a quantum well. 

Our paper is organized as follows: in section II we calculate the transmission and reflec- 
tion probabilities by using the transfer matrix technique; in section III we derive a formula 
for the density of states in the case of completely coherent transport. A completely anal- 
ogous formula which takes into account the effects of dissipation is obtained in section IV, 
and is shown to reduce down to the Breit-Wigner formulas for lower rates of incoherent 
processes in sec. V. A summary ends the paper. 



II. TRANSMISSION AND REFLECTION PROBABILITIES FOR DOUBLE 

BARRIERS 

Let us refer to the case sketched in Fig. 1: the one-dimensional potential energy profile 
V{x) defines the first barrier (a, 0), the well region (0, w), and the second barrier {w, b). Let 
us also introduce the wave vector k{x) = [2m{E — V{x))Y^'^ /h for all x where V{x) < E, 
where m is the electron effective mass in the material of the well, and h is the reduced 
Planck's constant. 

We can calculate the total transmission and reflection coefficients by using the transfer 
matrix technique.&^^ In the assumption of coherent transmission through each single bar- 



rier, the transfer matrix Mi for the first barrier satisfies all the properties listed in Ref. ^ 
and has the form 
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Ml 



(0 ^(0* /+W 



r?/t? 



(0* 



(1) 



where t^i and rf^ are the transmission and reflection coefficients, respectively, for a plane 
wave coming from the left electrode, with wave vector ku = k{a). The corresponding 
coefficients t^{^ and r^i^ for a particle coming from the right with wave vector kir = k{0) are 



t^i^ = ti'kir/kii and r\ 



(r) 



-r 



(I)*, (I) ,,(/)* 



t'-'/t'f'*. Moreover, if Ti = = \t'('>\'^kii/k 



^(012; 



Ir IS 



the transmission probability, and -Ri 



I'l I 



is the reflection probability, we have 



Ri + Ti = 1, i.e., the continuity equation for the probability density current holds true. The 
same considerations apply to the second barrier and its transfer matrix M2, provided that 
we define k2i = k{w), k2r = k{b) and change all the subscripts 1 into 2. 

In the well region, dissipative processes are accounted for by means of the mean free path 
/; that is, the intensity of a plane wave of wave vector k' has a decay length equal to /. As a 
consequence, the probability density current for a given wave function is not conserved. The 
effect of / is taken in to account by using a complex wave vector k[ = k' + i/2l: a plane wave 
of wave vector k[ along the x-axis has the form exp(i/c-a;) = exp(iA;'x) exp(— a;/2/), therefore 
its square modulus decays as exp(— x//). 

The multistep potential approximationil can be used to obtain the transfer matrix 
of the well, provided that the complex wave vector ki{x) = k{x) + i/2l is used at any point 
X in the well. If we make the hypothesis that theDotential varies smoothly enough that a 
semiclassical approximation is valid, we obtainciiiHB 
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where 
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exp t j ki{x)dx 



(2) 



(3) 



The transfer matrix Mdb for the whole double barrier structure is given b; 
Mdb = M,M^M2 



^/ '^db ' db /'' 



I ' db / ''db I'^db ' db' db / ''dbi J 



db 

mjr) /,(0l 



(4) 



where and r^^^ are the transmission and reflection coefficients, respectively, for an electron 
coming from the left, and t^^^ and r^^^ are the corresponding coefficients for an electron 
coming from the right electrode. Straightforward calculation yieldsil 



,(0 



.(0 

'^db 



''I '^2 7 
1-c 



and 



,(0 

db 



1 - c/Ri 



(5) 



(6) 



where 
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c = rl 7 . (7) 

The expressions for t^^"* and r^^"* can be easily obtained from (||) and (H) by substituting 
the subscripts / with r, 1 with 2, and viceversa. An electron coming from the left has a 
probability T^^ = |t^2P^2r/^ii of being transmitted and a probability -R^2 = kdfeP of being 
reflected; there is also a fraction S^^ = 1 — [T^j^ + r12) of electrons which have been absorbed 
(i.e., have undergone incoherent processes). In a completely analogous way we can define 
the corresponding probabilities T^l\ R^h, and S^f^ for an electron coming from the right. 
We also obtain that t]^'' = t]^-* — hence we will often write simply T^b — , while, in general, 



III. DENSITY OF STATES IN THE CASE OF COHERENT TRANSPORT 

In this section we will address the case of no incoherent process in the well, i.e., the mean 
free path I —>■ oo. In this case we have R^^l = R^^ = ^db = 1 — T^b and the results of Ref. |1| 
can be straightforwardly applied. 

We showedB that the density of states in a given system is equal to the sum of dwell 
times corresponding to each incoming channel divided by Planck's constant. In our case the 
region Vt of interest is {a,b), and there are two incoming channels, the left and the right 
ones, so that the density of states pn{E) in Q, including both spin contributions, can be 
written as 

Pn{E) = ^[4\E) + rj^\E)], (8) 
nil 

where r^^ and t^^^ are the dwell times for an electron of energy E coming from the left and 
the right electrode, respectively. 

In order to obtain rj^^ to substitute in (P), we can use the additivity of transmission 
and reflection times r^'' and t^^ obtained by using the Larmor clock and other well known 



approaches.liinl3 If we consider an electron coming from the left electrode, apply a uniform 
perturbative potential A on the double barrier {a < x < b), and re-calculate the total 
transmission and reflection coefficients as a function of A, we can write 



r«^Re{a|lu*S' 



Tjj'^ = Re I z^TT^ In 



dX 



(9) 

A=0 



(10) 

A=0 



and, finally, obtain0lli t^^ = r^Tdb + Rdb- In the Introduction we mentioned the 
controversy on the tunneling time problem, and we are aware of the fact that there is no 
wide consensus in the scientific community on the "actual" significance of Tj) and t'^ . 
Anyway, re-assuring the reader that we do not want to forget about the long debate in this 
field, for convenience reasons we will refer to (P) and (|10D as transmission and reflection 
times. 



5 



Substitution of (|) and (H) in (^0), after straightforward but cumbersome calculations, 
yields 

'''id = 'T'Tl^db + TRlRdb + TwTdb + 'Trl'^db 

+TrtT^-^^, (11) 

1 — cr 



where r^], Tj^g^ are the transmission times for barriers 1 and 2, and the well, respec- 

tively, defined as rjf-* in replacing with ti\ i li respectively; r]j'] is the reflection 
time for barrier 1, defined as in (|ToD with rf"* in the place of r^2- We call Trt round trip 
time] it is defined as 



r,, = Re I I = 41 + 2r^ + r^; (12) 

the last equality derives from (|^ and explains the name given to Trt'- it is actually the sum 
of the times corresponding to the steps needed for a round trip of the well: reflection from 
barrier 1, traversal of the well, reflection from barrier 2, and again traversal of the well. We 
can easily obtain r^"* by repeating all the passages from (|D to (pTD commuting the subscripts 
/ with r, and 1 with 2. If we substitute (|ll]) and the corresponding result for into (H) 
we obtain 

Pn{E) = \[{4lTab + ^Rdb) 
TxTi 

+ {4lTdb + 4lR 



db) 



db 



C 
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+ r.t^^(Ti+T2)]; (13) 

Pn(E) includes all the states in the region {a < x < b). We are actually interested in the 
states in the well region and in the tail states penetrating both barriers on the well side, i.e., 
the states in the "effective" well region, therefore we drop from pQ the terms which take into 
account the states on the left side of barrier 1 and on the right side of barrier 2 (i.e., the 
first and the second term of {^^, respectively). Moreover, the third term is easily shown to 
be much smaller (under the condition Ti, T2 <^ 1) than the fourth one, therefore the density 
of states Pw in the effective well region can be written as 



1 Ic 



2 



p^{E) = —Trt-^{T, + T2). (14) 
TTIl i — C 



From d^) and from the fact that Ti,T2 ^ 1 we have |cp(Ti — T2) ^ 1 — |cp, so we get 

1 



Pw = ^rrt-F(c), (15) 



where we have defined 

1 - Ic 



F{c) - ^ej^. (16) 



The density of states in the effective well region is therefore shown to be proportional to the 
round trip time times a factor F(c), which will be shown in the next section to depend only 
upon transmission and reflection probabilities for the whole structure and for each barrier. 
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IV. DENSITY OF STATES IN THE PRESENCE OF INCOHERENT PROCESSES 



A. Local density of states in the well 

In this section incoherent processes are taken into account, therefore the formula (P) for 
the density of states is no longer apphcable. The total density of states p^, in the effective well 
region is obtained as the sum of the density of tail states pi\E) and p2\E) penetrating 
both barriers on the well sides (Sec. IVB) and the integral of the local density of states 
p(x, E) in the well. 

In this section, in particular, we obtain a formula for p(x, E) which does not require 
the hypothesis of smooth potential in the well. We consider a point x' inside the well 
(x' G (0,w)). Let us split the x-axis into two regions, and let us consider the potentials 
Vi{x) = V{x) for X < x' and V{x') otherwise, and Vr{x) = V{x) for x > x' and V{x') 
otherwise, as sketched in Fig. 2. Let us call ri{x') the reflection coefficient for a plane wave 
of energy E incident on Vi{x) from the right, and rr{x') the reflection coefficient for a plane 
wave of energy E impinging on Vr{x) from the left. The local density of states at a point x' 
can easily be written asa 

where both spin contributions have been considered, the wave functions are not normalized 
and Jn{E) is the total current associated to state ipn entering the whole system. The sum 
is over all degenerate states corresponding to the same energy E, i.e., in our case, the ones 
associated to a particle coming from the left electrode (n = 1) and a particle coming from 
the right electrode {n = 2). The quantities to be put in (pIT]) are derived in the Appendix. 
Substitution of (|A1[ - |A^) and the corresponding quantities for ip2 in (P^ yields: 

TTh [ 1 — rr{x jri{x j J 

An identical result has been obtained through different procedures and in simplified condi- 
tions by other authors.i'S 



B. Density of tail states penetrating both barriers on the well sides 

For obtaining the density of tail states penetrating both barriers on the well sides we 
just have to use (|1^, provided that the sum is only over the states incident on the well side 

(r) 

of the barriers. Let us consider the first barrier: the density of tail states pi is 

We can find for p^^^ a more compact expression: in the Appendix we derived the currents 
j[ll^{x, E) and ^2mc(^' ^) incident on Vi{x), associated to the states ifji and il)2, respectively. 
Therefore, if we remember that the dwell time is defined as the ratio of the integral of the 
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probability density over the considered region to the incident current, from ( |A3| - [A4|) we can 
write 



nil 



(r) 



.h{E) 



J2{E) 



r, 



(r) 
Dl 



^rSFK(0)r,(0)], 
Tin 



(20) 



where the function F has been aheady defined in (|1( 

For the second barrier, following the same procedure, we have 



pI\e) = —T'j;>,F[rr{w)ri{w)] 



1 



(0 



(21) 



C. Density of states in the effective well region 

The density of states in the effective well region is therefore, from (|18D , (pOD, and (|2l|) , 

pUE) = p'f^E) + pI\e) + r p{x,E)dx. (22) 

We can write p{x, E) in a different way, in order to derive a more compact formula for 
Pw{E). Let the density matrix g{E) in the well be the incoherent superposition of states ipi 
and ip2 with probabilities pi and p2, i.e., 

{x\g{E)\x)=p,\Mx)\'+P2\Mx)\'. (23) 

Associated to g{E) there is the probability density current J{E), whose expression is given 
by ( [AID the corresponding quantity for ip2, that can be split into a left going com- 
ponent Ji{x,E) = pij[''^j^^{x, E) + p2J2inc{x, E) and a right going component Jr{x,E) = 
Pij[ll^{x, E) + p2J2il.c{^, E). Now, we can make the hypothesis that both Ji and Jr are 
much greater than the net current J = Ji — Jr- By imposing Jr ~ Ji we can obtain Jr, Ji, 
and {x\g{E)\x) as a function oipi/p2. This result, substituted in (p!8|), yields 



A great simplification of (|2^) can be obtained if we make again the semiclassical approx- 
imation in the well. In fact, in this case we have just to notice that rr{x')ri{x') = c for all 
x e (0,to). Therefore we can write 

pUE) = ^K,F{c) (25) 



where we have defined 



rrt-rni+rn2 + ^l J^x') + J.(x') ' ^ ^ 
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r^j can be interpreted as the round trip time of the well in the presence of inelastic processes. 



We wish to point out that (^5|) is formally analogous to (^) found in the case of coherent 
transport. It is also easy to verify that when / oo, i.e., when the limit of coherent 
transport is approached, r*^ tends to the value of defined in (p!2[). 

The condition Ti, T2 -C 1 allows us to write, with very good approximation, the following 
expression for F{c), where Tab, Rdb , and R^jl are obtained from (|) and (P): 



R 



Tdb 



R 



(r) 
db 



To 



Tdb 



(27) 



therefore 



can be written as 



nil 



rC. 



(0 

■db 



Ti 



+ 



T 



db 



(28) 



clS db function of the round trip time and transmission and reflection probabilities for the 
total structure and for each barrier. 



V. COMPARISON WITH BREIT-WIGNER FORMULAS 

In this section we want to show that for a rate of incoherent processes low enough (i.e., 
a long enough mean free path), the formula for the density of states derived above reduces 
to the one obtained by the means of Breit-Wigner formulas. 

Let us expand c given by (^ to first order around the resonant energy (which is the 
energy at which c is real and positive): 



c{E) ^ c{Er) 



'^—^{E-En) 



(29) 



where we have used the definition (^) of and the fact that a:Tg{ihd log c/dE} ^ 0. If 
w/l is small enough we can write 

c{En) = [R,{En)R2{En)Y/'e--/' 

^ 1 - ^T,iEn) - l^iEn) - j. (30) 



Substitution of (|2|) and (|0|) into (|) and (D yields 



-db 



(31) 



{E - Enf + (r/2)2 
and 

1 _ ^(0 _ ri(r2 + rj) , . 

^ ^^^^ [E-EnY + {V/2r ^ ^ 

where Ti = hT,{En)/4,{En), T2 = hT2{ER)/4,{ER), = h2w/[4,(ER)l], and T = Ti + 

r2 + r,. 
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Equations (|3TD and ( p2D are the Breit-Wigner formulasjilii and Fi, F2, and Fj are 
the partial resonance widths for each process allowing escape from the resonant state, in 
particular tunneling through barrier 1 and 2, and incoherent processes, respectively. Partial 
resonant widths are characteristic quantities of the Breit-Wigner formalism, and are given 
by the ratio between % and the characteristic time of the process we are considering. In 
the case of escape through one of the barriers the characteristic time is intuitively given by 
the ratio of the round trip time and the tunneling probability of the barrier. In the case 
of inelastic scattering the time is r^t times the ratio between the mean free path I and the 
length corresponding to a round trip of the well (2w). 

From (^8]) and (|3T| - |32D , we straightforwardly have 

1 F 

P""^^^ '^{E-ERy + {Vl2r ^^^^ 

i.e., the result usually obtained from Breit-Wigner formulas.iB We wish to point out that 
this formula holds true if the development of c to first order oi E — Er and to first order 
in w// is a good approximation. In other words, Breit-Wigner formulas can be used if each 
partial width is much smaller than both the resonant energy E^ and the difference between 
the height of the barriers and E^. In our case these conditions are true for Fi and F2 and 
holds true for F2 if //w is high enough. Otherwise, the expression given by (p5|), which has 
a wider range of applicability, has to be used. 



VI. SUMMARY 

In this paper we have studied the density of states in a double barrier structure. We 
have proposed a simple model which is able to account for inelastic processes occuring in 
the quantum well by means of a single phenomenological parameter /, the mean free path. 

We have obtained a very compact formula which relates the density of states in the effec- 
tive well region to the round trip time of the quantum well and to the tunneling probabilities 
for the single barriers and for the whole structure. The formula is valid both in the case 
of completely coherent transport and in the case when dissipative processes in the well are 
predominant. 

This formula will be shown to be fundamental in unifying two widely known description 
of transport in double barrier structures: that of resonant tunneling and the one of sequential 
tunneling .0 

We believe that the role of the density of states — a characteristic quantity of a system 
in equilibrium — in the steady state transport and in the characteristic times of the motion 
of electrons in a mesoscopic system deserves a deeper investigation. 
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APPENDIX: 



Let us consider a point a; = a;' in the well and the state ijji of energy E corresponding 
to a particle coming from the left of x'. We can describe ipi as a plane wave of amplitude 1 
undergoing multiple reflections on Vr{x) for x > x' and Vi{x) for x < x' (see Figs. 1 and 2), 
so that we have 

V^l(x') = 1 + Trix') + Trix'Mx') + ■■■= , ^ ^ J'jf] ■ (Al) 

1 - rr[x')ri[x') 

Now, [1 — |r;(a;')p] is the probability that a particle impinging on Vi{x) is not reflected back 
(i.e., is either transmitted or "absorbed" on the left of x'). For time reversal symmetry, it is 
also the probability that an electron coming from the left of x' appears at x — x': if 1 is the 
amplitude of ipi before taking into account multiple reflections, the total current entering 
the system has to be 

1 - \ri{x')\^ 

where v{x) = hk{x)/m. It is worthy noticing that the dependence of Ji on x' is due only to 
the fact that ipi is not normalized. We can also associate to ipi and x' a probability current 
density which can be split into a left going component Junc (incident on Vi(x)), and a right 
going component jillc (incident on Vr{x)): 

. /^ _ \rr{x')\Mx') 

-^^inci^)- \l_r,(x')rJx'W ^ ^ 



and 



4ncM- \l-ri{x')lr{x')\^' ^"^^^ 

The corresponding quantities for the state 1^2, associated to a particle coming from the right 
oi X — x', can be obtained by substitution of 1 with 2, r with /, and viceversa. 
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FIGURES 
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FIG. 1. The one-dimensional potential V{x) defines the two barriers [(a, 0) and (u),5)] and the 
well (0,iy). 
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FIG. 2. The double barrier structure is split into two regions at x = x': rr{x') is the reflection 
coefficient for a plane wave of energy E incident on Vr{x) from the left and ri{x') is the reflection 
coeffincient for a plane wave of energy E incident on Vi{x) from the right. 
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